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O 1 ABSTRACT 

O ' 



In dilute astrophysical plasmas, the collisional mean free path of a particle can 
exceed its Larmor radius. Under these conditions, the thermal conductivity and 
viscosity of the plasma can be dramatically altered. This alteration allows out- 
wardly decreasing angular velocity or temperature gradients to become strongly 
destabilizing. This paper generalizes an earlier, simple analysis of the viscosity 



instability, by including the dynamical effects of magnetic field line tension. Such 
\q ', effects lower the growth rates found in the absence of such tension, but still al- 

low growth rates in excess of the maximum of the standard magnet orot at ional 
instability. We find very good quantitative agreement with more complex kinetic 



treatments of the same process. The combination of large growth rates and large 
magnetic Prandtl number suggest that protogalactic disks are powerful dynamos. 



O 

6 . 
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1. Introduction 

Magnetic fields, even when highly subthermal, turn free energy gradients in sufficiently 
ionized fluids into sources of instability. There are important astrophysical consequences 
of this result: by the magnetorotational instability (MRI; Balbus & Hawley 1991), ac- 
cretion disks become turbulent when angular velocity Q, rather than angular momentum 
QR 2 , decreases outward, and by the magnetothermal instability (MTI; Balbus 2001), dilute, 
stratified plasmas are destablized when temperature rather than entropy decreases upwards. 
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Since it is easier to violate the free energy gradient criteria than the classical Rayleigh and 
Schwarzschild criteria, it is the former, when applicable, that are relevant to the behavior of 
their host systems. 

The MRI is the best-known example of a free-energy gradient instability, and has by now 
been extensively studied. Recently, however, it has been noted that purely viscous effects in 
a dilute MHD fluid can also lead to accretion disk turbulence, even if the J X B Lorentz 
term (J is the current density and B is the magnetic density) is negligible (Balbus 2004; see 
Quataert et al. 2002 and Sharma et al. 2003 for plasma kinetic treatments). This is because 
there is a large parameter domain in which the ion Larmor radius is very small compared 
with fluid length scales (bounding the field strength from below), yet the magnetic field 
contribution to fluid stresses is tiny (bounding the field strength from above). If the plasma 
in question is sufficiently dilute, the ion cyclotron period is much smaller than the ion-ion 
collision time. Under these circumstances, the viscous stress tensor is highly anisotropic, and 
dominated by its diagonal terms (6). These circumstances are similar to those of anisotropic 
thermal conduction. In the MTI, heat, borne by the electrons, flows only along the magnetic 
lines of force. In the viscous case, it is of course angular momentum flow that is restricted. In 
both cases, fluid attributes that are ordinarily responsible for stabilizing dissipation (thermal 
conduction, viscosity) become active agents of destabilization, a remarkable turn of events. 

The cause of the instability is essentially the same in both cases. Perturbed magnetic 
field lines, initially isothermal (isorotational), are stretched along the direction of the back- 
ground temperature (angular velocity) gradient. This allows heat (angular momentum) to 
flow from one fluid element to another. The elements then move yet farther apart, the field 
lines become more aligned with the background gradient, and the preocess runs away. The 
magnetothermal instability has been followed calculated far into the nonlinear regime, in 
which vigorous convection has developed (Parrish & Stone 2005). 

Neglect of the fluid J x B force in the magnetoviscous calculation of Balbus (2004), 
while enabling a point of principle to be made more transparently, ultimately restricts the 
domain of validity. Magnetic stresses are an important complication, acting as both a sta- 
bilizing and destabilizing agent. To understand the true nature of this instability in an 
astrophysical context, it is desirable to pursue a more general approach, and this is the 
primary motivation of the current paper. 

In §2, we discuss the physical parameter regime in which this work is expected to 
apply, and present the formalism by which the components of the magnetic viscosity may be 
calculated. In §3, we formulate the problem and obtain the dispersion relation of interest. 
Characteristic growth rates are determined, and the physical regime reexamined for self- 
consistency. In §4, we present a discussion relating our work to a more complex plasma kinetic 
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treatment, and touch upon some astrophysical implications. Finally, in §V, we summarize 
our conclusions. 



2. Preliminaries 
2.1. Parameters 

The magnetoviscous instability finds its natural venue in protogalactic disks and halos, 
as well as in very low density accretion flows, of which the Galactic center is the prototype 
(Melia & Falcke 2000). These systems are characterized by subthermal magnetic fields in 
collisionless plasmas and ohmic diffusion coefficients small relative to the viscous diffusion 
coefficient. Let u C i be the ion cyclotron frequency, and Tj the ion-ion collision time: 

- - 6 - 91 x 105 (sex) ^ 8 - w 

where n is the proton density in cm -3 , T4 the ion kinetic temperature in units of 10 4 K and 
In A is Coulomb logarithm. We shall work in the asymptotic domain 

1.09 x 10 5 ^ T 3 4 /2 B,g ^ x 



n J In A 

where B^ G is the magnetic field strength in microgauss. With n < 1 and T 4 > 1, even very 
weak fields can be accommodated by this regime. On the other hand, we shall also assume 
that the Reynolds number is large, which requires 

ftn < 1, (3) 

where Q is the disk rotational frequency. The collision frequency 1/Tj is therefore much 
larger than an orbital frequency, but much smaller than the cyclotron frequency. Finally, 
let us note the Prandtl number ratio of the viscosity to ohmic resistivity (Balbus & Hawley 
1998): 

We shall assume V >> 1. The Spitzer (1962) viscosity for a hydrogenic plasma is: 

77 ~ 1.1 x 10- 16 T 5 / 2 (J^ g cm- 1 s" 1 . (5) 

(Applications to low luminosity black hole accretion should use a value of In A closer to 30.) 
Finally, for future reference, we use the standard plasma (3 parameter to represent the ratio 
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of the gas to magnetic pressures: 

? = (•) 

It is helpful to have representative physical parameters at hand, even if they are very 
crude. For protogalactic disks, densities could range from 10~ 2 to 1 particle per cm 3 ; tem- 
peratures from 10 4 to 10 6 K. For radiatively inefficient black hole accretion flows, a typi- 
cal density might be 10 8 cm -3 , but the range of interest should be thought of as perhaps 
10 5 — 10 10 . The ion and electron temperatures seem to be very different in these flows: the 
ions ought to be virialized at T ~ 10 12 K, whereas the electrons are thought to be ^ 10 10 
K (Narayan, Mahadevan, & Quataert 1998 for a review). In a forthcoming paper, we shall 
discuss a possible reason for this two temperature structure based on the MTI. 



2.2. Magnetic Viscosity 

A convenient formalism for the viscous stress tensor in the presence of a magnetic field 
is presented in Balbus (2004), and shall quote the results here for reference, referring the 
reader to this paper for further details. 

Our fundamental coordinate system for the disk will be a standard cylindrical system: 
radius R, azimuth 0, and axial variable Z. Define now a local Cartesian system, determined 
by the magnetic field. We denote these axes by subscript b. Z b points along the local 
direction of the magnetic field, X b and Y b may be any axes orthogonal to one another as 
well as to the Z b direction. Following Braginskii (1965), the Z b Z b component of the viscous 
stress <Jz b z b is unaffected by the presence of the field. It is given by (Balbus 2004): 

cr Zb z b = -277 [(6 • V)v] -b, (7) 

where v the local velocity field, and b is a unit vector in the direction of the magnetic field. 
The other diagonal components of the traceless viscous stress are 

°x b x b = ov b Y b = ~(T ZbZb /2. (8) 

To find the components of the magnetized viscous stress tensor in any other local Cartesian 
frame, the transformation law may be written 

a a = U-jb) °i b j b , (9) 

ibJb 

where once again the b subscript denotes the magnetic field frame and bold face quantities are 
unit vectors of the indicated component. Using equation (8) for the nonvanishing diagonal 
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stress tensor components in the field frame, we obtain 



- cr Zb z b 



(i-Z b )(j.Z b ) - l -{i-Y b ){j-Y b ) - l -{i.X b ){j.X b ) 



(10) 



Once <Tz b z b is determined in the magnetic field frame, may be calculated in any frame. 



3. Formulation of the Problem 

3.1. Braginskii Stress 

We consider the stability of a disk under the influence of a weak magnetic field. As in 
standard MRI analyses, we assume that the field is sufficiently weak that it has no effect on 
the equilibrium state, but that magnetic tension is important for the behavior of local WKB 
perturbations. The fundamental fluid equations used here are mass conservation, 

H + v.( pv ) = o, (11) 

the equation of motion, 

" (I + "' v ) • = - v ( p + 8 + tJ b ■ V > B - " v * - % < 12 > 

and the induction equation of ideal MHD, 

— = Vx(vxB). (13) 

$ represents an external gravitational potential and the other symbols have their usual 
meanings; the viscous stress tensor is given by equation (10). 

The equilibrium state is a differentially rotating disk. As stated above, we work in a 
standard cylindrical coordinate system, R, 0, Z. The angular velocity is f2(i?), and we shall 
restrict ourselves to a local analysis at the midplane. Thus, we may ignore buoyant forces. 

In the equilibrium state, it is assumed that = 0; it will be shown that this is 
in general not a stable configuration. The initial magnetic field lines are spooled around 
cylinders, unaffected by the shear. 

We consider next small departures from the equilibrium flow. Linearly perturbed quan- 
tities are denoted by dv,5aij, etc. We work in the local WKB limit, with the space-time 
dependence of all perturbed quantities given by exp^t + ik-r). Thus, 7 is a growth or decay 
rate if it is real, and an angular frequency if it is imaginary. 
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To evaluate Sa^ we follow Balbus (2004). Since vanishes in the equilibrium state, it 
follows from equation (10) that 



i-Z b )(j-Z b ) - l -(i.Y b )(j-Y b ) - l -(i-X b ){j-X b ) 



Sa i:j = 5a Zb z b 

The geometry of the equilbrium field is defined by 

b = cos x </> + sin x Z. 



(14) 



(15) 



4>, Z, and R (used below) are unit vectors in the indicated cylindrical directions, and x 
is the angle between the magnetic field and the axis. The local magnetic field axes are 
chosen to be 

Z b = b = cos x 4> + sin x Z, (16) 
X b = RxZ b = -sinx + cosx Z, (17) 

and 

Y b = R. (18) 

See fig. 2. 

The diagonal elements of Sa^ are found to be 



5<?rr = --S(rz b z b , 



cos x- 



sm x 



S°z b z b , fozz = sin 2 x - 



cos 2 x 



5(J z b z b , 



and the off-diagonal elements are 

3 

Satz = Sa z<t> = -cosx sinx Sa Zb z b - 
Finally, we evaulate 5az b z b from its vector- invariant form (7), 

6a Zb z b = ~2v ([(Sb-V)v] -b + [(b-V)Sv] -b + [(b-V)u] -Sb) 
Using v = RQ<fi and equation (16), we obtain 

dtt 



5a 



z b z b 



-2r] 



5b 



R 



d\nR 



cosx + i(k-b)(b-Sv) 



(19) 
(20) 

(21) 
(22) 
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3.2. Linearized Equations 



The linearized dynamical equations are 

k-dv = 0, 



ik 



-fSv R - 2Q6V.P + —[5P- 



5a 



ZbZb 



p \ 2 

2 3K 



+ 



B ■ SB 



4:71 



ik z B z sinx 
4np 



SB R = 0, 



2fi ^^/T Sm X C ° S X ^ Attp = 

(Here, k 2 is the square of the epicyclic frequency, k 2 = 4f2 2 + df^/dlni?.) 



SP+ I sin 2 0-^) <)a 



ZbZb 



ik z B z sin x 
4np 



= 0. 



The induction equations are 



^5B R - ik z B z sin x&v R = 0, 



7 5B* - 5B R 



k z B z sinx^v = 0, 



d\nR 

j5B z — ik z B z sin x$ v z = 0. 

Finally, we recast equation (22) for the all-important Z b Z b stress component in 
convenient form, 



8az b z b = -2i]ik z sin x 



cos x | SV ^ R + Svj, ) + sin x <foz 



(23) 
(24) 

(25) 
(26) 

(27) 

(28) 
(29) 
a more 

(30) 



3.3. Dispersion Relation 



The linearized equations lead, after straightforward but somewhat tedious algebra, to a 
dispersion relation that blends the properties of the magnetorotational and magentoviscous 
instabilities: 

k 2 ( k 2 \ 

^ 7 4 + 3r] vl 3 sin 2 X k 2 ± (k 2 R + k z cos 2 X ) + 7 2 + 2^- (k • v A f j + 



(k • v A ) k ± + k z cos x 



dfl 2 ' 
dlnR 



k 



da 



k 2 z {k ' VA)2+ dlnR 



0, 



(31) 
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where k\ = k\ + k\ cos 2 x is the square of the component of the wavenumber perpendicular 
to b, and i] V is the kinematic viscosity i]/p. We note that in the limit r\y — > the standard 
MRI dispersion relation is recovered, while in the limit k • — > 0, the magnetoviscous 
dispersion relation of Balbus (2004) emerges (taking also the limit k^ = 0). Unstable modes 
are present if dQ 2 /dR < 0, as expected. 

To evaluate the growth rates associated with the dynamical magnetoviscous instability, 
we proceed as follows. First, we set k^ = 0, which corresponds to the most rapidly growing 
modes, and simplifies the calculation. Second, we calculate all rates (7, k • v^, Vvk±) in 
units of Q. We introduce the notation 

X = (k • v A ) 2 /n\ Y = 3r] V sin 2 x^lM 7 = 7M (32) 

and restrict ourselves to the astrophysically interesting case of dfl 2 /dR < 0. Our dispersion 
relation may then be rewritten as 

7 /4 + 7' 2 [(t/^) 2 + 2X] + X(X + d\nn 2 /d\nR) 
' " Y{\dhiQ 2 /dlnR\ -X-y 2 ) ' ^ > 

Triplets of the form (X, Y, 7') are easily calculated by this formula, and contour plots of 
7 thereby generated. In fact, it is convenient to view the results in the Y/X, X plane 
since the wavenumber is thereby removed from the ordinate, which becomes a normalized 
viscosity. The results are shown in figures 1 (Keplerian rotation profile) and 2 (galactic 
rotation profile). The maximum possible growth rate is 



7 2 

1 max 



dtt 2 



d\nR 



(34) 



as found by plasma kinetic treatments (Quataert et al. 2002; Sharma et al. 2003) and for 
the magnetoviscous instability (Balbus 2004). The effect of magnetic tension is to lower the 
maximum growth rate to somewhere between this value, and the MRI maximum rate of 
0.5|dn/dlni2|. 



3.4. Physical Regime 

Figures (1) and (2) show that magnetoviscous effects can significantly change the inviscid 
MRI growth rates in the parameter regime X < 1, F > 1. Let us determine what sort of 
magnetic field strengths are implied by this condition in an environment representative of a 
protogalaxy. With sin 2 x set equal to a half, 

Y = = ±^W) 2 , (35) 
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Keplerian Rotation 




0.5 1 1.5 2 2.5 3 

(k-v A /Q) 2 



Fig. 1. — Contours of normalized growth rate 7/fi for a Keplerian rotation profile, as a func- 
tion of normalized wavenumber (k • v^/fl) 2 and normalized viscosity 3^f2cos 2 x/ v a- Con- 
tours run from 7 = 0.5f2 (contour a) to 7 = 1.5Q (contour b) in steps of 0.1Q. The maximum 
MRI growth rate is 0.75fi. 



where we have introduced the isothermal sound speed eg and scale height H defined by 
HQ = cs- Replacing rjv by pi] and using equation (5), we find 

y= 2.39x10-1^ 

n 
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Galactic Rotation 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 

(k-v A /Q) 2 

Fig. 2. — Contours of normalized growth rate j/tt for a galactic rotation profile (Q oc R^ 1 ), 
as a function of normalized wavenumber (k • v^/fi) 2 and normalized viscosity 3r)Q cos 2 x/ v a- 
Contours run from 7 = OAQ (contour a) to 7 = 1.2Q (contour b) in steps of 0.08O. The 
maximum MRI growth rate is 0.5f2. 

where T 4 is the temperature in units of 10 4 K, and Q the angular rotation rate in units of 
10~ 14 rad s _1 . In other words, 

nT' 3/2 

(k ± H) 2 £ 4.2 x 10 7 — 4 , (37) 

''-14 



for viscous effects to be important. On the other hand, X < 1 for the growth rates to be 
significantly enhanced. Since 



X = 




(k z H) 




(38) 



we find that 



P > 2k 2 z H 2 > 8.4 x f 0' 



7 




) 



(39) 



(ignoring the distinction between kz and k±.) This translates to a magnetic field strength 



This is of the order of the strength of the seed field estimated in Balbus (2004), obtained by 
diluting stellar surface mangetic fields of strength ~ 0.1 G over parsec scales. For applications 
to inefficiently radiating black hole accretion flows, the permissible f3 parameter regime can 
extend to values near unity. 



Anisotropic stresses in the presence of a magnetic field may be calculated from a moment 
expansion of the collisionless Boltzmann equation using a formalism developed by Kulsrud 
(1983). (In the absence of a heat flux, this reduces to a double-adiabatic equation of state 
first derived by Chew et al. [1956].) This approach has been applied to the MRI in the fully 
collisionless regime (Quataert et al. 2002) and the collisional regime (Sharma et al. 2003; 
see also Snyder, Hammett & Dorland 1997). In the fully collisional limit, Sharma et al. 
recover the MRI, as do we in the limit rjv — > 0. The fluid and plasma kinetic treatments 
are markedly different in their level of complexity, and a quantitative comparison of their 
findings is therefore of interest. This is shown in figures 3 and 4. We have adopted the 
fiducial case of Sharma et al. (2003). This corresponds to x — 45°, (5z = SnP/B 2 , = 10 4 . 
Choosing values for k ■ and v/Vt, where v is the ion self-collision frequency, then allows 

our X and Y parameters to be determined, once use is made of the Spitzer viscosity (5). 
The figures show the growth rates as a function of k ■ v A /il for several different values 
of v/VL. In general the two approaches are in excellent agreement (growth rates agree to 
within ~ 5), with significant disparity coming only in the completely collisionless regime 
(infinite viscosity). This is not unexpected, since it corresponds to the complete breakdown 
of the fluid regime. What is remarkable is how well the relatively simple fluid treatment 
performs. The savings of effort is the greater the more complex the problem. For example, 



of 



B & 6.5 x 10- 10 T 4 5/4 Q^ 



G. 



(40) 



4. Comparison with Plasma Kinetic Treatment 
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Fig. 3. — Dispersion relation for the magnetoviscous instability based on a kinetic treatment 
(figure adapted from Sharma et al. 2003). 



an investigation combining both viscous and thermal effects is relatively straightforward via 
a fluid treatment; both effects need to be included in a study of high temperature accretion 
flows. This will be presented in a forthcoming paper by the authors. 

It should be noted, however, that the effort entailed in a kinetic treatment has valuable 
dividends. It is difficult to know, for example, how far into the long mean free path regime 
an MHD fluid approach may be relied upon before collective instabilities might be triggered. 
This level of analysis can be provided only by a kinetic treatment (e.g., Schekochihin et al. 
2005). 
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1.8 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

k z v Az /Q 

Fig. 4. — Dispersion relation for the magnetoviscous instability using an MHD fluid ap- 
proach. 

5. Conclusions 

We have analyzed the stability of a dilute astrophysical plasma in the presence of Bragin- 
skii (1965) viscosity. Our main result is the dispersion relation (31). The regime of interest 
corresponds to gases in which the ion gyroradius is small compared to a mean free path. In 
laboratories this means a very strong field, but in astrophysical applications even very weak 
fields are in this regime when the density is sufficiently small. This regime does not apply 
to classical accretion disks or to stars (the densities are far too high), but it does apply to 
the interstellar medium of galaxies and protogalaxies, and to the low density accretion flows 
around some giant black holes (e.g. the Galactic center). 
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The work presented here extends the study of Balbus (2004), which included Braginskii 
viscosity but ignored J X B Lorentz forces, and complements the plasma kinetic treatment 
of Sharma et al. (2003). The dynamical effects of the magnetic field lower the maximum 
growth rate found in the absence of such effects (eq. [34]), but growth rates significantly in 
excess of the MRI value 0.5\dQ/d\nR\ are still found. 

The magnetoviscous instability thus enhances the dynamo effect of the MRI in dilute 
gases, and therefore may well be a significant source of magnetic field amplifcation. Whether 
it is a more efficient mechanism than others that have been suggested (e.g. Schekochihin 
et al. 2004) is yet to be established, but numerical simulation of this large Prandtl number 
regime is clearly now a viable possibility. At the very least, MRI growth rates are significantly 
enhanced in the linear regime, and nonlinear field dissipation is inhibited (since the resistive 
scale is "blocked" in a downward cascade by dissipation at the larger viscous scale). 

In the presence of strong temperature gradients in accretion flows, magnetothermal and 
magnetoviscous effects may both be important. Since a mere temperature gradient can now 
trigger convective instability, consequences for the temperature structure of inefficiently ra- 
diating accretion flows can in principle be profound. These will be discussed in a forthcoming 
paper. 
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